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We consider the influence of optical and temperature-dependent atomic fluctuations on the forma- 
tion and propagation of optical vortex solitons in dense media realized as hollow-core optical fibers 
filled with a cold atomic gas in presence of optical pumping. We show different perturbation-induced 
scenaria of complete destruction and smooth transformations of the topological characteristics of 
localized optical structures in hollow-core fiber. The maximum levels of optical and atomic fluctua- 
tions at which the soliton regime can be maintained has been determined. The estimates for these 
levels show an opportunity to observe the optical vortex solitions in the core-filling gas of the fiber 
for temperatures smaller than the critical temperature for Bose-Einstein condensate. 
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I. INTRODUCTION 

The study of the formation and propagation of local- 
ized spatial optical structures in different media 1] is 
one of the important issues of modern applied atom op- 
tics. It is associated with a wide range of applications, 
some of which are the problems of data transmission and 
processing Q. The greatest interest here is the special 
class of optical topological structures known as optical 
vortices [3|. The dark center for those could be regis- 
tered reliably in experiments even in the case of a strong 
diffraction broadening of the optical beams It is rel- 
atively easy to obtain such structures in laser cavities 
@ or spatially inhomogeneous media @. However, the 
maintenance of the stability conditions for such struc- 
tures when they propagate in the real media is a difhcult 
task. The soliton modes of optical vortices have so far 
been observed quite rarely even within short distances 
0. At the same time the distances required for telecom- 
munications have not been worked through at all. The 
explanation for this is the fact that a full-scale solution 
to the problem of developing such information channels 
requires a correct determination of the stability regions 
Q fulfilled in the parameters of a laser experiment as well 
as the stress tests with perturbations for the propagating 
localized optical structures. 

Formation of the stable optical structures occurs near 
the lasing threshold [2], where the nonlinear effects are 
very strong. Therefore the solution of this problem re- 
quires a strict account of the higher-order nonlinearities 
and nonlinear absorption effects which are crucial for sta- 
bility of the optical structures 2, 9] . In experiments in 
order to realize an optical control of the parameters when 
entering the soliton mode it is most beneficial to use the 
A-type interaction configuration in a resonant atomic 
medium with a probe localized structure and a pump 
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field. The observations of the strong nonlinear effects in 
such a scheme can be associated with the realization of 
the Raman regime [1, [l3| , when the detunings of the elec- 
tromagnetic fields from the atomic transitions are much 
higher than the corresponding decay rates (with account 
for the optical depth of the medium, see [ll|b 

Among promising media for observation and control of 
solitons we can mark out the recently created hollow-core 
optical fibres filled with a cold atomic gas [l3|. First, it 
would allow us to observe the coherent Raman scattering 
and study the competition between the nonlinear, diffrac- 
tion and dissipative effects for the probe pulse over short 
distances [H, [l3| . Second, the use of cold atoms opens 
new possibilities of precision and efficient control over 
the optical properties of the system. Indeed, in a cold 
medium we may neglect the effects of the spectral line 
broadening and splitting as well as assume the minimum 
impact from the thermal noise. 

Obviously, the quality of transmission and the process- 
ing rate in the optical data channel jisj is strongly de- 
pendent on the efficiency of nonlinear transformations, 
which significantly increases with the density of reso- 
nant atoms in the medium. However, for media with 
a large density of optically active particles it is necessary 
to take into account the effective or the local value of 
the field acting on them The atom-field interac- 

tion becomes trickier in this case and in the literature it is 
known as the consequence of the near dipole-dipole inter- 
actions (NDD). This may significantly modify the picture 
of the phenomenon as the conventional Rabi frequency 
for an atomic transition is comparable to the correspond- 
ing NDD correction term. The estimations for ensembles 
of resonant particles show that the NDD or the local 
field effects may be very significant starting from den- 
sity p = 10^^ cm"'^. In connection to this complication, 
some special attention should be paid to the problem of 
obtaining stable dynamics of spatial solitons and optical 
vortices in particular [3] in dense media with an addi- 
tional view to the perturbations in the medium and the 
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field. 

The influence of the noise effects on soliton solutions of 
the nonlinear Schrodinger equation (NLS) has been well 
studied in the literature and a violation of dynamic equi- 
librium in such a system is usually observed for the values 
of the noise intensity, comparable to the main parameters 
of the problem |23| ■ A particularly strong presence of the 
noise in this problem manifests itself near the zero dis- 
persion point, when the short-term changes of the disper- 
sion sign initiated by fluctuations in the medium would 
quickly destroy the soliton. Only some original pinning 
method can help to restore the system's working char- 
acteristics [21'] . Analysis of the fluctuations influence on 
the dissipative solitons is usually performed artificially 
by introducing a noise source to the Ginzburg-Landau 
equation (GLE) [2^ [2^. However, this approach does 
not answer the question about the nature of the noise, 
and is therefore not possible to assess the contribution of 
the various physical components into the overall picture 
of the development of instabilities in the system. This 
is especially true in the case of the optical solitons in 
gaseous media, where the density fluctuations at small 
scales can be comparable with the average density over 
the entire medium [l^, ■ Even transition to the Bose- 
Einstein condensate (BEG) can not solve the problem 
completely because in such a system there are intricate 
dependencies of density fluctuations on the trapping ge- 
ometry and the type of interaction between the particles 

Mm- 

In this paper we consider the problem of obtaining 
a stable spatial dynamics of optical vortices Q in the 
dense media of BEG-filled optical fibers in presence of 
optical pumping. The problem is similar to the observa- 
tion of temporal self-induced transparency solitons in mi- 
crostructured materials considering the influence of the 
local field effects Secondly, in this paper we study 

the effect of random and/or periodic perturbation of the 
system's parameters on the process of stabilization and 
shape evolution of the vortex solitons propagating in a 
fiber filled with a cold gas. 



II. GENERALIZED DENSITY MATRIX 
EQUATIONS FOR A-TYPE ATOM-FIELD 
INTERACTION IN A DENSE MEDIUM 

In this paper we assume that a probe light pulse of 
given shape Ep with the center frequency propagates 
along axis z of a hollow-core optical fiber filled with a 
gas of cold atoms with a temperature T near the critical 
temperature [2^ of the phase transition: 



where h is the Planck constant, fcs is the Boltzmann con- 
stant, N is the total number of particles, ujq is the trap- 
ping frequency, and C(3) is the zeta function of Riemann. 
It propagates in the direction opposite to cw pump radi- 
ation Ec as illustrated in Fig.lTJi. It propagates in the di- 
rection opposite to cw pump radiation Ec as illustrated in 
Fig. [1^. In the Raman limit the probe field detuning 
should be substantially greater than the relaxation rate 
of Tab i^ac) from the excited state |a) (see Fig. [T)d). As 
applied to this problem, the parameter of optical depth 
of the medium do can be determined through the char- 
acteristic linear dimension oq of topological structures 
formed in the xOy plane as do — 9^Nao/{cTac) (com- 
pare with ^10]). Here g = ^ba y^u!/2heoV is the atom- 
field coupling constant, e = Ap^/2eoV/hw, Ap is slowly 
varying amplitude of the probe field, V is the quanti- 
zation volume, N = pVo is the number of atoms in the 
interaction zone of volume Vq for atomic density p, c is 
the speed of light in vacuum, Sq is the vacuum permit- 
tivity. The frequency separation between levels \b) and 
|c) is (5 = 6.834 x 10^ s~^, the dipole transition matrix 
element for \b) \a) is pba = 3.58 x 10"^^ G-m. 

Since ^^Rb atoms have large dipole moments of the 
desired transitions and taking into account their density 
inside the hollow core of the fiber the value of NDD in- 
teraction Xba — P \Pba\^ /(3feo) which determines the ef- 
fective Rabi frequency i^eff = ^ba + Xba^ba [3 '^^^ be 
close to the value of the Rabi frequency of the probe field 
(here aba is the corresponding of the density matrix ele- 
ment) . Thus, the nonlinear optical effects arising in such 
a system will be determined, on one hand, by the pump 
field parameters, and on the other hand, by the number 
density of the atoms hosted by the hollow fiber. 

This system is described by the following system of 
equation for the density matrix elements: 




Figure 1. (a) Model diagram of a hollow-core optical fiber 
filled with an atomic gas (b) A-type interaction configuration 
applied for ®^Rb atoms, energy level splitting between |c) and 
j6) is (5 = 6.834 X 10^ s~^, the dipole transition matrix element 
for \a) \b) at wavelength A = 780.241 nm is ^ab = 3.58 x 
10"^^ C-m. 
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where VL = fica and ge = ^Iba are the Rabi frequencies 
for the pump field and the probe pulse respectively, At 
and Ac are the corresponding detunings. 

In this section, we consider the applicability of ap- 
proximations used in the paper for a three-level atomic 
medium as a result of the different modes of its interac- 
tion with the external optical fields. One of the specific 
requirements to such a problem in our case is the pres- 
ence of a significant polarization of the medium produced 
by the contiguous transitions \a) — )■ \b) and |a) — t- |c). It 
is necessary to provide nonlinear control over the probe 
vortices via the pump wave. Besides, we have to make a 
set of assumptions that will significantly simplify Eqs. ([T]) 
for further analysis. The key point here is to choose the 
right balance between the Rabi frequency of the pump 
field, its detuning from the resonance and the decay rate 
of the excited state. 

In general, the propagation equation for a probe field 
in a resonant medium, taking diffraction into account, is 
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where V^^ — d'^ /dx^ + d"^ jdy^ , D = A/tt the diffraction 
parameter in the plane transverse to axis z. 

Self-consistent solutions of Eqs. ([ll)-® can be found 
if we use the following approximations. First, we as- 
sume that all atoms are initially at level 5, i.e., abb = 1, 
faa — <^cc — and the population of the excited state re- 
mains small throughout the interaction time, i.e., Ubb = 
1, (Jaa(cc) — (while &ii = where i = a,b,c). Secondly, 
we believe that when Xca = Xba = X the contribution 
of the local field is comparable to the probe Rabi fre- 
quency, but it can be omitted for the pump transition, 
i.e., ge > (JbaXba and Q. ^ cFcaXca- As a result, the final 
system of equations reduces to the following (simplified) 
form: 

o-ba = -riCTba - ige - iflcTbc - iXba<Jba, (3a) 

CTca = -^2<yca - iS^^CTcb " iXbaOba^cb, (3b) 

dbc ^ i^ac^bc + igscTac - i^*crba, (3c) 
fJbb = igeaab - ig£*o'ba, (3d) 

where Ti = -iAb-Hl/2(r„b+rac), T2 = -zAc-hl/2(r„b+ 
Tac) and A3 = Ab — Ac. In Eqs. ^ the contribution 
of the local field is presented in two ways. In ([5a|l its 
action is trivial and produces an effective frequency shift. 
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Figure 2. Time dependence of the density matrix elements in 
units of the excited state lifetime F";,^ (a) for the near-resonant 
(b) Raman modes of A-type interaction in the limit of strong 
coupling, (c) for the near-resonant mode in the limit of weak 
coupling. The solid line is for \aba\^, dashed is for \abb\'^ and 
the dash-dot line is for |crccp, the red lines are solutions of ((T}, 
the dotted blue line is for |af,a|^ obtained from Eqs. (|3}. The 
thick solid green line corresponds to |(Ti,aP which is calculated 
for the approximate steady-state solution ((§]). 



but in (|3bp it provides a significant effect, which is the 
appearance of the nonlinear coupling between the atomic 
excitations (polarizations) in the probe aba and magnetic 
(Tcb transitions. 

Figure [2] shows, in particular, the comparison of the 
independent solutions of the exact system ([l} and the ap- 
proximate ^ for different modes of the atom-field inter- 
actions at probe field switching time of Tsw ~ 2x 10^^" s 
and cw pump field. The atomic and field parameters were 
chosen as follows: the characteristic size of the optical 
beam ao = 20 ^m, the atomic density p = 1.01 x 10^^ 
m~^, the relaxation rate F^b = Fac = 10^ s~^, the probe 
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field intensity Ip = 0.22 W-cm~^, the pump intensity 
Ic = 146.5 W-cm^^, the probe detuning Af, = — 5 x 10^ 
s^^. The corresponding Rabi frequencies can be calcu- 
lated as = iiacEc/h and ge = fJ.abEp/h using the field 
strengths -Ec(p) — \/Ic{p)'^/ceo and be = 1.13 x 10^" 



s~^ and ge = 4.4 x 10® 



respectively (g = 1.3 x 10° 



s~^ and N = 27.7 x 10®). These parameters meet the 
condition Tab < ^ for the strong atom-field coupling in 
the system. 

At near-resonance interaction, when Ac/doTac ~ 0.1 
at do = 0.308 and Ac = 3 x 10^ s~^, the nonlinear cou- 
pling between the probe and pump fields is actually ab- 
sent as it follows from Fig. [51 This is due to the fact 
that during the interaction the ground state \b) is hardly 
depopulated {abb ~ 1) and in the absence of the atoms in 
states I a) and |c) the resulting polarization conditioned 
by transition \a) — \b) is weak (see in Fig. [5^). In this 
case, the solutions of Eqs. Q and (jS)) for the matrix el- 
ement aba can approximate one another with sufficient 
accuracy as a proof of validity of approximations used in 
the derivation of ([3]). Figure [2^ shows the regime typical 
for the linear FIT effect in a three-level medium [2^, or 
its non-linear analogue [30j at a significant increase of the 
pump intensity. 

Quite different dynamics related to the nonlinear inter- 
action of the fields may be expected when we chose the 
Raman interaction mode when Ac > $7 > doTab- This 
interaction mode is long-known [Toj and well studied, 
but it did not get a strong "scientific resonance" in terms 
of its prospects of practical use in modern optical tech- 
nologies, as it happened with FIT. This mode is char- 
acterized by a set of curves for the matrix elements in 
Fig. obtained for parameters the same as in Fig. [^^ 
but with greater detuning Ac = 3 x 10"'^° and condi- 
tion Ac/doTac ~ 100. In this case the atomic state |6) is 
effectively depopulated giving population mostly to state 
|c). This presupposes the emergence of a significant po- 
larization in transition |a) — )• It seems apparent if 
we compare the scale of values a^a in Figs. [5^ and [^h 
(the same is true for the transition \a) — |c)). This fact 
leads to the establishment of a nonlinear energy transfer 
between modes of the probe and pump fields. This is 
the interaction mode we are going to use further in our 
analysis. 

The weak coupling in the atom-field system Tab > ^ 
for both modes in Figs. [2^ and [2]d produces a signif- 
icant discrepancy between solutions obtained from ([T]) 
and ^ as one may see in Fig. [2t at Ac = 3 x 10^ s~^, 

= 1.13 X 10® s~^. In particular, this is due to a viola- 
tion of the condition ft >• acaXca and the need to account 
for the effects of the local field for transition \a) — >■ |c). 
This fact complicates the task and requires a direct nu- 



merical simulation of the self-consistent problem ([2]) and 
([3]) [m . It makes further analysis of opportunities to ob- 
tain solitons, including spatial solitons, very troublesome 
because of major transformations of the matrix elements. 

Now let us find the steady state solution for the density 
matrix element aba of the probe transition in a form that 
depends only on the material parameters of the media 
and the characteristics of the optical fields considering 
various cases of atom-field interaction. In order to do 
this, we should solve Eqs. ^ in two steps. The first 
step is to define the polarization of the system at the 
lower levels acb as based on the approximation of constant 
populations and atomic polarizations in the steady state, 

i.e., CTaa = lytb = (Jcc ^ and aba = <^ca = O-fec = 0. At 

this point from Eqs. ^ we obtain an algebraic equation 
for the polarization abc- 



T,g'\er 



the roots of which are the following: 



CTfcc = 



r;A) abc + ig£^*T*2 = 0, (4) 



(5) 



where D = - {r^g^ \e\' + TlAf ~ Aixg^ Tl\e\\ 

A^\n\^-iA^{T^+tx), x = Xba- 

Solutions ^ determine, in fact, the two branches of 
spin excitations that occur at the transition between lev- 
els 1 6) and |c) (see Fig.[T]). The solution that contains the 
minus sign in ([5|) leads to a problem with the saturating 
nonlinearity abc ~ 1/e and is not considered in this pa- 
per. Expansion of the other solution, that is positive, in 
series of the pump field s gives the following relationship: 
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Our second step is to use the equations for both a^a and 
abc in (IS]) to write the atomic polarization from the probe 
transition: 
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Now we substitute expansion (jS]) into ([7]) and perform 
the secondary expansion to get 
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(8) 



In Fig. [5] solution is represented by a thick solid line, 
which is in agreement with the solutions of ([ij and ([3]) 
to a high degree of accuracy as a clear proof of the cor- 
rectness of ([8]). Let us note that if solving Eqs. ([3]) in 
the EIT limit, when it can be assumed that aac = 0, 
the expression for ate contains only the first term in the 
right-hand side of (|6]) (see [295). its turn, Eq. ([7]) now 
becomes 



AaCTbc 5A3 
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(9) 



and determines the appearance of the phase modulation 
and the absorption (amplification) of the probe pulse in 
the atomic medium. We found that the appropriate solu- 
tion for this linear mode Q for aba with sufficient accu- 
racy coincides with the solutions of ([T]) and in Fig. [5^. 
The later proves the validity of the theory and the limit- 
ing cases as well. 

For further analysis we will focus on the situation 
shown in Fig. [^b and move on to the study of the spa- 
tial dynamics (in the plane perpendicular to the direction 
of propagation) of optical pulses in the hollow-core fiber 
filled with resonant atoms as shown in Fig. [U and con- 
sider the possibility of obtaining a special kind of stable 
structures known as vortex solitons 3] of the probe field. 



III. VARIATIONAL APPROACH FOR STABLE 
VORTEX SOLITONS IN A THREE-LEVEL 
MEDIUM 

In the Raman limit for the A-type interaction scheme 
(Fig. [1]) after substituting ([8]) into the propagation equa- 
tion ([2]), the self-consistent problem of spatial dynam- 
ics (H])-® is reduced to the well-known form of the 
Ginzburg-Landau equation (see Q): 
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-ai£ — a2\£\ e — a4|£| e. 



where the corresponding coefficients are the following: 
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is the cubic nonlinearity 

^ f ig^N I 1 
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is the quintic nonlinearity; 

ai — Im 

is the linear loss coefficient 



a2 — Rc ■ 



Ac 



^1 ix , -xM^ 



\ A A* 



A^ 



is the cubic loss; 
ai = Re< — 



(ATlY 



A3/^2 ^X'\^ 



^X 




is the quintic loss. 

It is obvious that Eq. (fTO|) is the result of a compli- 
cated nonlinear interaction between the probe field and 
the pump field which arises solely due to the presence 
of significant values of polarizations of the optical tran- 
sitions in the Raman limit of the A-configuration jTlj . 
Even though the locally acting field used in this the- 
ory does not produce new terms in Eq. (jlOp it, however, 
gives certain corrections to the picture of the nonlinear 
interaction (terms with x~fa'Ctor) for the case of an op- 
tically dense medium. Let us note that in another limit 
of the A-configuration the near-resonance conditions, i.e. 
Af, <^ Tab, would simplify Eq. ([TOl) by eliminating its 
nonlinear terms. The local field in this case would be 
insignificant as it can only lead to an additional phase 
modulation in the probe field 32 1. 
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Figure 3. Frequency dependences of the third 72 (red solid 
hne) and fifth-order 74 (blue dotted line) nonlinearities, and 
the backward diffraction (green dashed line) . The inset shows 
the enlarged image of the region near the point A. The inter- 
action parameters are the same as in Fig. [2jb). 



concept of dissipative solitons, to maintain the energy of 
bright sohtons at a constant level it would require inter- 
changing between absorption and gain effects in different 
areas of the probe envelope, in particular, when these in- 
equalities (5>0, (/)<0, /i>0 are to be met [34|. Figure[4] 
shows the typical dissipative parameters as functions of 
the pump field intensity Ic at a fixed detuning as 
in Fig. [3] near the gain threshold. Under conditions of 
diverse variation of the dissipative coefficients with in- 
creasing pump intensity we can select Ic for fixed (or 
vice versa) to satisfy the condition of dynamic equilib- 
rium for nonlinear and dissipative processes. 

Using the variational approach for the analysis of (fTTj) 
we may determine the ranges of parameters for which 
stable spatial solitons can occur. Our particular interest 
is focused on an important class of dissipative spatial 
vortex solitons of the kind described in : 
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In order to start the analysis of ((TOl) we will first trans- 
form to the moving coordinate system T = t — z/c 

and perform the change of variables u = e / \J \ein^ , 

^ = z/Ldf, X = x/qq, Y ~ y/aQ. Our second step 
is to define the following main characteristic lengths: 
■^72 = 1/(72 kinl^) and L^4 = 1/(74 |£m|^) are the 
nonlinearities of the third and the fifth orders respec- 
tively; Lai — 1/a/ describes the linear losses, La2 = 
l/(a2|em|^) and L^a = l/(ci!4 |£in|^) describe the non- 
linear losses of the third and fifth orders; Ldf — a^/D 
is the diffraction length, where is the reduced ampli- 
tude of the beam at the entrance point. After multiplying 
both sides of pO)) by Ldf we finally get 



dU 



dY^ J 



\u\^u- 



-,.\ufu^Q, (11) 



where Q = i[—6U — (j)\U\^ U — jx |i7|** t/] is the dissipative 
part, for which we have the following notations U = uN 
and N"^ — Ldf/L^2, and the main characteristic param- 
eters: S = Ldf /Lai, L^2/La2, M = ^^2/(^04^^/), 

= ^72/(^74^^/)- 

For the formation of a conservative optical soliton 
the interaction parameters should be selected in such 
a way that the characteristic nonlinear and diffraction 
lengths were approximately comparable or multiples of 
each other [s^. This can be achieved near the region 
of strong transformations of the nonlinear coefficients by 
varying the detuning of the pump field as it is shown in 
Fig. [3] In particular, for the point A^ the diffraction 
length should be only Ldf = 1.61 mm so for the positive 
nonlinearity, i.e. 72 > 0, we may expect focusing and for- 
mation of a stable spatial soliton. However, following the 



where r — Va^ + F^, 9 is the angle in spherical coordi- 
nates and A, R, C, and ^ are the amplitude, the spatial 
width, the curvature of the wave front and the phase 
of the pulse, respectively. Parameter S determines the 
topological charge of the vortex soliton. In the special 
case of 5 = 1 the normalization coefficients Aq and Rq 
can be expressed in terms of the full power normalization 



P = 



\U{r,9)\^rdrd0 = (ttSIAIrI) A^R^ 



In the simplest case, when P = A R , they are related 
as Aq = l/{Roy^) and Rq = 1. 

In practice, the effects of diffusion in atomic ensembles 
would fill the dark center of an optical vortex with light 
([T^ even with presence of a long-lasting shape in the 
form of a bright soliton propagating in such media [3]. 
For a self-maintained special crater-shaped vortex soli- 
tons (jl2p the system would require not only a balance of 
the nonlinear, diffraction and dissipation effects 3J| but 
also some presence of the optical diffusion 35| or a com- 
plicated type of modulation of either the refractive index 
|36| or the absorption coefficient Q or both. In partic- 
ular, the static parameter of the linear absorption can 
be replaced with a space-dependent effective parameter 
<5cff = S — Vr^ . It would describe introduction of an addi- 
tional saturable absorber into the fiber filled with a cold 
gas. This is the case we consider further in this paper and 
assume V = —0.03 everywhere below. A detailed analy- 
sis of the stability of vortex solitons in the linearization of 
the master equation written in a form more general than 
([T0| and with the optical diffusion taken into account is 
given in 



37| 



Now we use the Euler-Lagrange equation for the ax- 
isymmetric vortex soliton and we get the system of equa- 
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Figure 4. Parameters of the dissipative part of the non- 
dimensional Ginzburg-Landau equation 5 (solid red), (p (dot- 
ted blue), and ^ (dashed green) as functions of the pump 
field intensity. The interaction parameters are the same as 
for point A in Fig. [S] 



tions for the variable parameters: 
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dR _ (j)A^R 2fiA'^R 
dC" ~ IGtt ^ 8l7r2 ^ 
1 



(13a) 

CR + VR^, (13b) 
2iyA^ 



dC ' 8i?4 167ri?2 ' 8l7r2i?2' 

_ 3^2 iQj^^i I 

2R^' 



dC 8tt 



8l7r2 



(13c) 
(13d) 



For a low frequency modulation {C^ ^ 0), as in [sj, we 
come to the following system of equations, which simpli- 
fies the search for the fixed points for 



^2 (-8l7r(/) - 32A2^\ 



12967r2 



^2 (SlTT - 32^2^,) ' 



1627r2 



A2 (8l7r - 32^2^,) ' 



(14a) 
(14b) 



167r 8l7r2 



+ A{-6-C + VR'^) =0. (14c) 



Equations ([T^ give 16 steady state roots, among which 
only two correspond to the physical conditions on the 
energy and the width of the vortex soliton (A > 0, i? > 
and A,Rg 3?). Besides, only one of them would have the 
maximum value of A and negative frequency modulation 
(C < 0) and be stable Q. 

Figure [5] shows the parametric plane formed by the 
following parameters: the density of resonant atoms in 
the system p and the intensity of the pump field Ic- The 
region bounded by the dash-dot line refers to stability of 



axisymmetric vortex solitons and arise for the selected 
physical solution of ([T^ Q. This area was determined 
from the analysis of the eigenvalues of the Jacobi matrix 
for Eqs. ([T^ . i.e., for condition Re{Xj) < 0, where j — 
1, 2, 3 (38|], and corresponds to the point of a stable focus. 

Direct numerical simulations of Eq. pT|) with variation 
of parameters for the function of the form ([T^ and taking 
into account the initial angular perturbation for R and 
C [39] show that the area of stable solutions obtained 
from simplified Eqs. (1141) is very approximate. At the 
same time the true stability region I (in Fig. [5]) is much 
smaller. The region of stability obtained using the vari- 
ation approximation appears to have a "fine" structure 
in the form of separate zones for solitons with modified 
shapes as well as zone IV, where the optical vortices are 
damped. 

In particular, in the area II in Fig. [5] (parameters as 
in Eq. v = 0.1653, 5 = 0.6819, (j) = -1.1998, 

^l = 0.2432 at p = 9.85 x 10^^ m'^ and h = 146.2 
W-cm~2 for point B) we find a spontaneous transition of 
an axisymmetric vortex into a single-humped asymmet- 
ric stable vortex soliton, presented in image B in Fig. |6l 
Bifurcations of this type are described in Q. Region 

III in Fig. [S] (parameters as in Eq. dTT]): v = 0.2519, 
5 = 0.7158, (j) = -1.2576, ^l = 0.3169 at p = 9.85 x 102i 
m^"^, Ic = 152 W-cm~2 for point C) is characterized by 
a high sensitivity of the system to initial perturbations. 
As these perturbations are developing the vortex soliton 
spontaneously loses its topological charge and transforms 
into a new stable state with S = (see image C in Fig.[6|). 

In area V the account for angular perturbation in 
Eq. p2|) leads to destruction of vortex solitons and emer- 
gence of nonstationary localized structures in their place, 
which, however, do not damp but demonstrate some on- 
going evolution [40]. At point D the field parameters in 
Eq. ^ take values v = 0.1043, 6 = 0.674, </) = -1.14, 

= 0.1551 at p = 9.6 x 102^ m"^ and Ic = 143.3 
W-cm~2. Here the axisymmetric vortex soliton trans- 
forms in a multi-hump structure rotating along its propa- 
gation path in the fiber as it is shown in image D in Fig. [51 
It is also possible to observe a specific double-hump soli- 
ton js^ in image D' in Fig. [51 This structure occurs 
when v = 0.1297, 6 = 0.6633, (/) = -1.1848, ^ = 0.228 at 
p = 9.925 X 1021 m-3 and Ic = 143.3 W-cni-2. In area 

IV we notice the loss of stability and complete damping 
of vortex solitons. 

Let us note that point A (with the same set of param- 
eters as in Fig. is indeed found in the area of stability, 
calculated using both the variation approach and direct 
simulations, which is in agreement with the qualitative 
analysis of Eq. PT|) . Thus, formation of vortex solitons 
in our medium would require the initial spatial width and 
the curvature of the wave front of the light pulse in the 
form of Eq. ^her^^ 19 pm and Cr = -3.7x10^ m-2 
respectively. The dimensionless parameters are, there- 
fore, A = 5.75, R = 0.95, and C = -0.0015. The 
values calculated for Eq. ([Tl]) for point A in Fig. [5] are 




9 9.4 9.8 10.2 

px10'\m-' 

Figure 5. Parametric plane (the intensity of the pump field 7c 
vs the dopant density p). The dash-dot hne margins the area 
of stability for an axisymmetric vortex soliton obtained by 
the variation technique. The highlighted areas are obtained 
by direct numerical simulations and correspond to: green - 
axisymmetric vortex solitons, yellow - single-peak vortex soli- 
tons, red - irrotational solitons, orange - unsteady localized 
structures, purple - unstable structures. The interaction pa- 
rameters correspond to Fig. 2] 



V = 0.1899, 5 = 0.675, </> = -1.2375, and /i = 0.3075. 
The soliton regime of such a structure is achieved vifithin 
the characteristic length oi Lst = 4 cm while the ro- 
tation period of the vortex wave front in this regime is 
Tv — 1.45 cm (see image A in Fig. [B]). 

Figure [7] shows the region of stability for vortex solitons 
in the plane (/c, V) a.t p = 1.015 x 10^^ m~^. The area of 
stability obtained by the variational approach extends to 
both the positive and negative areas of parameter V. The 
numerical simulation gives a narrower area and demon- 
strates that stable dissipative vortices can exist only for 
negative values of V. 

The fundamental point in this problem is considera- 
tion of the local field effects. Indeed, if x = 0, that is, 
without taking them into account, the region of stability 
is not simply transformed, but it completely transcends 
the parametric areas shown in Fig. |4l In this case, all 
the solutions obtained for the vortex solitons become un- 
stable. Thus, the strategy of the possible experiment 
should change dramatically with the local value of the 
acting field in the medium. 

Further analysis and evaluation of the other stabil- 
ity regions for the vortex solitons require solution of the 
full nonlinear system (jl3p as well as an expanded multi- 
dimensional numerical experiment within the domain of 
the system's parameters. 
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Figure 6. Results of direct numerical simulations of Eq. 
i.e. spatial profiles in {X, F)-plane of optical beams at the en- 
trance ^ = and after traveling the distance of ^ = 100000 in 
the gas-filled fiber in the presence of azimuthal perturbations. 
The letter in the upper left corner of each image corresponds 
to a point in the parameter plane in Fig. (5] the coordinates of 
which were used to calculate the parameters of the Eq. (jlip . 

IV. STABILITY STRESS TESTING OF 
OPTICAL VORTICES UNDER PERTURBATIONS 
OF SYSTEM'S BASIC PARAMETERS 

Let us consider the evolution of an optical vortex soli- 
ton (fT^ in terms of natural fluctuations or/and addi- 
tional modulations of the optical fields and the atomic 
medium in the core of the fiber. Such perturbations may 
be associated with variations of the pump field intensity 
Ic{^,X,Y) = i;?"* + Q{^,X,Y) and the atomic density 
p{£.,X,Y) = + Cp(C,X,r), where I^"^ and p'^''* de- 
scribe the deterministic part. Parameters C/(Ci Y) and 
(^p{£,, X,Y) are presented by three-dimensional matrices 
of fluctuations or periodic spatial modulation of the cor- 
responding functions. Random perturbations of the in- 
tensity Q'^'^{^,X,Y) we describe as the Gaussian white 
noise. The density fluctuations in a Bose gas below the 
critical temperature at transformation to the dimensional 
parameters have a spatial correlation function of the form 

(Cp"'('^i)Cr'(^2)) - Po{S{n - r,) + i.„(r,„d)), (15) 
where 

In our numerical simulations we used the three char- 
acteristic spatial scales, i.e., the sampling interval of 



154 



150 



O 



146 



142 
-0.10 




(a) /„Wcm-^,. 




Figure 7. Parametric plane (tiie intensity of the pump field 
Ic vs. the optical trapping parameter V^) . Interaction param- 
eters and filling colours are the same as in Fig. [S] 



^st = 0.0011, the characteristic length of density vari- 
ation (physical step) of ^rnd — 0.0068 (equivalent to 
''rnd = 11 /ini) and fvar = 0.07 for the characteris- 
tic distance at which the mean value of the noise be- 
comes zero. The selected physical step is so great that 
'^ni'i^i-nd) / Po ~ 0.4%, so the second term in the right-hand 
side of Eq. is neglected and we can model the den- 
sity fluctuations using the white noise and spatial splines 
(Fig. IH^jb) for the convenience of calculations at the in- 
termediate points. 

The periodic perturbations of the pump field inten- 
sity CT^iO the atomic density Cp^^iO '^ith the 
amplitudes Cil^ and Cpo^ are assumed to be effective 

only along axis z: Ci'^iO — C^^ ^^^^(2^?/^/ + Vi) and 
Cp^^iS,) = sin(27rC/Lp + ipp) with the specified spatial 
periods Lj, Lp and random initial phases fi{p) ■ The pul- 
sations of the pump intensity can have two sources: the 
artificial modulation at the output of the pump laser or 
the intensity oscillations occuring in the pump depletion 
and gain saturation f4l| . The periodic modulation of the 
density can be caused by an external mechanical action 
entailing excitation of acoustic waves in the system. 

The specificity of our problem lies in the possibility of 
studying the overall picture of the development of per- 
turbations of each GLE parameter and their impact on 
the dynamics of the soliton based on the perturbations of 
the basic parameters describing the media and the con- 
trol fields. Quantitatively, this may be reflected by the 
values of the noise strength Dn = Zvar(^fi)^ for the GLE 
coefficients (72, 74, a/, 02, 04) with their mean-square de- 
viations, calculated based on the presence of random per- 
turbations in the system £7 = cr///^^* and Sp = ap/p'^°^, 



Figure 8. The model of the statistical stability for the axisym- 
metric vortex solitons with fiuctuations of the atomic density 
p (a), the intensity of the pump field Ic (b). The range of 
the fluctuating quantities between the red lines corresponds 
to the stability area in Fig. [5] The temporal dynamics of the 
reduced power I of the soliton is shown in image (c), the pro- 
files of the wave packet in image (d). The letter at the upper 
left corner of each image corresponds to the time in image (c) . 



where ct/ and <Tp are the mean-square deviations of vari- 
ables, and Zvar = £,va.i- Ldf- The possibility in our numer- 
ical simulation of collective impacts of noise -I- periodic 
modulation of the basic parameters completes the picture 
of the full-scale stress testing in the propagation regime 
of spatially localized structures. 

To study the effect of perturbations on the vortex dy- 
namics we performed a series of numerical experiments 
for different values of aj and cTp and the control param- 
eters 1^°*^ and p'^°*^. As an indirect criterium of stability 
we assume the statistical stability of the reduced power 



\U{r,e)\^rdrde 



shown in Figs. ISllTOl under perturbations X,Y) and 
Cp{S,, X, y), and then we restore and analyze numerically 
the shape of the soliton in the desired point on axis ^. 

The longitudinal perturbations give rise to the inert 
properties of the soliton, i.e., the presence of noise, in 
which the phase trajectories p(^, X, Y),Ic{£,, X, Y) never 
leave the stability region I in Fig. [Sj does not cause dra- 
matic consequences. Beyond the stability area the dy- 
namics of a vortex is determined by the retardation ef- 
fects in changing its shape compared to the rapid changes 
in GLE parameters initiated by the noise. Figure [8] 
shows the case of the maximum-allowable noise param- 
eters ef = 5.1% and e'^J' = 7.4% {I^"^ = 146.5 W-cm'^ 
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Figure 9. The model of the "breathing" state as a result of 
migration of the soliton solutions between regions III and V 
in Fig. [5] in presence of fluctuations (a) the atomic density 
matrix p and (b) the pump intensity Ic- Image (c) shows 
the temporal dynamics of 7; corresponding to different time 
profiles of the probe wave packet in image (d). 



and p'^'^' = 10.1 X lO'^^ m^'^ correspond to point A in 
Fig. [5]) for the vortex soliton which preserves its shape 
over long distances z = 16.1 m at ^ = 10000 for values 
of the noise strength Z?^ = 3.2 x 10"^, D^, = 1.1 x 10"^, 
Ds = 1.66 X 10"^, = 1.1 X IQ-^, see the inset in 
Fig. m The presence of the transverse perturbations in 
plane xOy leads to the noise-induced distortion of the 
wave packet profile in Fig. (B}: but has almost no effect on 
its stability. 

The continuous noise parameters in Eq. ()11|) for a point 
Jdet = W-cm-2 and p^'^* = 9.77 x lO^i m^^ from 
the transition area II in Fig. [S] can initiate reduction of 
topology and a jump transition of vortex solitons in a 
"breathing" mode with constant interconversion between 
the fundamental soliton and oscillating polygonal struc- 
ture as in Fig. [HI This evolution is the result of walks 
the optical structure performs between areas III and V 
in Fig. [S] for values of noise strength — 1.8 x 10^^, 
= 7.23 X 10-^ Ds ^ 1.2 x 10-^ = 7.4 x 10"^ 
which correspond to the perturbations £/ = 1.92% and 
Ep = 2.87% at ^rnd = 0.05. In transition of a soliton 
to the vortex-free mode one may observe formation of 
crescent-shaped optical structure as shown in Fig. I^Jl. In 
further rising of ct/ and Cp there comes a transition into 
an unstable mode and decay, which is associated with a 
disproportionate increase in the noise intensity for vari- 
ous GLE coefficients. In this model the fifth-order non- 
linear absorption coefficient p has a more rapid increase 




Figure 10. Model as a result of the loss of stability (a,b) 
periodic modulation of control parameters in the presence of 
a low-intensity noise. Image (c) shows the temporal dynamics 
of /; corresponding to different time profiles of the probe wave 
packet in image (d). 



in the noise level and brings the soliton to zone IV in 
Fig. [5] and destroys it there. 

The destructive effect are also caused by a weak peri- 
odic modulation of control parameters in the presence of 
a low intensity noise in Fig. [101 which collectively lead to 
fast violations of dynamic equilibrium in the system and 
damping of the soliton at scales of the order of ^ = 8. The 
simulation parameters correspond to = 7.5 W/cm^, 
ipi = -2.51 and Cp^^ = 0.75 x lO^^ m'^, ipp = 2.76, 
Lj = 6, Lp — 2 at constant values of /^°* and p'^°' as 
for point A in Fig. [5l The scale of density fluctuations in 
Fig. [TO] is 3.2 mm that given the expression for the speed 
of sound in a Bose gas Cgn = (/i/mo) V47raap [12] cor- 
responds to the frequency 1 Hz of the sound wave with 
its speed 3.2 mm/s, where toq is the mass of a single 
atom and Oa is the scattering amplitude. In the pro- 
cess of damping of the vortex soliton the system shows 
the rapid development of angular perturbations leading 
to the emergence of two or multihump unstable localized 
structures. This behavior is associated with a large spa- 
tial scale of the simulated pulsing perturbations which 
keep the vortex in violation of the conditions of dynamic 
equilibrium for a long time. 

The well-known estimates for the modern laser systems 
that are compatible with the telecommunication channels 
determine the maximum allowable noise level of the out- 
put intensity to be not higher than 1%. This estimation 
by a large margin falls within the considered theory of 
stability of vortex solitons. Particle number fluctuations 
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of a Bose gas in the range of T < Tcr have the tempera- 
ture dependence [i^ 



-^N fl- 
ea^) \Tcr 



This equation provides a simple expression for estimat- 
ing the maximum density fluctuation in BEC Sp ~ 



^(6C(3)7V) 



-1/2 



= 1.177V~^/^ under condition that T = 



T^r- For TV = 3 X lO'* and = 50 kHz the critical 
temperature is 72 /^K. The maximum density fluctuation 
in such a system is 0.68%, which is 10 times as low as 
the maximum allowed atomic noise calculated for soliton 
mode in Fig. [51 This fact makes if evident to observe the 
vortex soliton regime at a temperature below the critical 
temperature for BEC. 

In Figs. [51ITU1 we also consider the fluctuations of 
diffraction parameter emerging in the cross-section with 
the spatial scale within the range of the field's wave- 
length. This requirement arises from the fact that the 
present fluctuations of the gas density effectively modu- 
late the refractive index of the medium. They also create 
a random phase delay, which with the finite width of the 
emission line produce the wavelength noise parameter. 
However, the simulation shows a statistically stable pic- 
ture with the presence of local perturbations of soliton 
shapes even at 10% deviations of A. This is similar to the 
initial angular perturbations of the vortex shapes which 
are smoothed after elimination of the noise. 



V. CONCLUSION 

This study may give some ideas that could be use- 
ful for development of optical communication channels 
using vortex solitons in hollow-core fibers filled with a 
cold atomic gas. With the given fiber's parameters and 
the proper intensity of the probe field one may simul- 
taneously tune the intensity and frequency of cw pump 
field to provide the conditions for the soliton propagation 
regime of vortex structures in the optical system. Be- 
cause of the nonlinear contribution from the local field 
effects to development of competitive optical processes 
in a dense atomic medium the NDD corrections should 



adjust the values of the control parameters for which we 
expect appearance of the solitons. For the parameters 
used in direct simulations the maximum fiuctuations of 
the pump intensity and the probe wavelength that al- 
low the soliton regime impose minor restrictions on the 
stability of the pump laser power and the probe laser fre- 
quency. The estimates for the atomic fluctuations gives 
an opportunity to observe the optical vortex solitions in 
the core-filling gas of the fiber for temperatures smaller 
than the critical temperature for BEC. The large-scale 
perturbations in the development of acoustic waves and 
the modulations of the pump field in the fiber can quickly 
destroy the soliton propagation regime even if they are 
of small amplitude. 

For long-term maintenance of the BEC state in the 
information channel it is possible to use fine fibers filled 
with BEC atoms with a core diameter of the order of 
one micron or less |4^]. Longer lifetimes of a coherent 
state in such a system are due to the geometry of the 
fiber and, in addition, the effect of channeling of atoms 
at creating the surface light wave along its core, causing 
the atoms to lose their energy effectively. The estimates 
given in [Z^l show that for the density 10^^ cm~^ the 
temperature inside the fiber drops down to 1.5 x 10~^ K. 

However, such an elongated BEC must be assumed as a 
ID gas so it is necessary to use some different relations for 
the fluctuations [i^ at which observation of the vortex 
solitons is possible. Except it, the use of thin hollow-core 
flbers for optical communication due to the need to create 
localized optical structures at nanoscale, but in solving 
such a problem transition to the Maxwell equations for 
continuous media is impossible, and the effects of the 
near-field will play a crucial role. 
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